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We investigate bipartite entanglement in spin- 1/2 systems on a generic lattice. For states that are an equal 
superposition of elements of a group G of spin flips acting on the fully polarized state |0)® n , we find that the 
von Neumann entropy depends only on the boundary between the two subsystems A and B. These states are 
stabilized by the group G. A physical realization of such states is given by the ground state manifold of the 
Kitaev's model on a Riemann surface of genus g. For a square lattice, we find that the entropy of entanglement 
is bounded from above and below by functions linear in the perimeter of the subsystem A and is equal to the 
perimeter (up to an additive constant) when A is convex. The entropy of entanglement is shown to be related 
to the topological order of this model. Finally, we find that some of the ground states are absolutely entangled, 
i.e., no partition has zero entanglement. We also provide several examples for the square lattice. 
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I. INTRODUCTION 

Entanglement emerged recently as a quintessential concept 
in several fields of physics. In quantum information theory 
(QIT) entanglement is a sine qua non resource for various 
quantum processing and quantum communications protocols, 
like teleportation, dense coding, cryptography, and is crucial 
for the exponential speed-up of several quantum algorithms 
Q|. The same concept is also essential for our understand- 
ing of several solid-state systems. Examples are the entangled 
ground state for two highly non-classical systems: supercon- 
ductivity (the BCS state |2, 3]) and fractional quantum Hall 
effect (the Laughlin state |4]). Another manifestation of the 
ubiquity of entanglement is found in the study of quantum 
phase transitions, where it is believed to be responsible for the 
appearance of long-range correlations 01 . 

Spin systems are also a distinguished playground for the 
study of bipartite entanglement and its scaling with the sub- 
system size. It has been shown that systems in which entan- 
glement scales less than logarithmic can be efficiently simu- 
lated on a classical computer |6]. Hence the amount of entan- 
glement present in the system and its scaling is crucial for ef- 
ficient quantum algorithms, i.e., problems that are classically 
intractable. In the case of a critical spin chain in XY and 
Heisenberg models, the entanglement between a spin block of 
size L and the rest of the chain scales like S ~ log 2 L and thus 
this system can be simulated classically |7|. Several groups 
have analyzed recently the entanglement properties of various 
spin systems |8], including 1 -dimensional lattice models of 
the XY 001, Heisenberg (I and Anderson model iToll . 

In this article we study bipartite entanglement in general 
spin systems extending our previous results found for the 
ground state of the Kitaev's model 1 1 1]. The relevance of this 
model stems from the fact that it was the first example of the 
new subject of topological quantum computation I YA [1311 and 
because it features topological order This is a type of 
quantum order which describes states of matter that are not 
associated to symmetries, like the fractional quantum Hall liq- 
uids. We show that for a large class of states (e.g., states stabi- 
lized by groups of spin flips, which include the ground state of 



the Kitaev model) the von Neumann entropy depends only on 
the degrees of freedom belonging to the boundary of the two 
subsystems. Hence, our result echoes the holographic princi- 
ple 1151 : the geometric entropy of a region A depends only on 
the degrees of freedom of the boundary of A, and not of the 
bulk. 

The structure of the paper is the following. In Section II 
we expose the general formalism for bipartite entanglement 
in spin systems and we apply it to states that are stabilized 
by groups of spin flips. For spins on a lattice we provide a 
geometrical interpretation of these results. In Section III we 
exemplify this general framework for the Kitaev's model flUl 
and we apply it to calculate the ground state entanglement. We 
calculate analytically the entropy of entanglement for several 
partitions of the lattice, like spin chains and spin ladders in 
Section IV. We then conclude in Section V. 



II. ENTANGLEMENT IN A SPIN SYSTEM 
A. The reduced density matrix 

In this section we find a general expression for the reduced 
density matrix of an arbitrary spin system and we give a neces- 
sary and sufficient condition for its diagonality. These results 
specialize in a very interesting way to those states, like the Ki- 
taev's ground states (see SectionllTH that can be written 
as an equal superposition of all the elements in a group acting 
on the reference state |0). To begin with we do not assume 
any particular geometry or dimensionality of the spin system. 

Given a system of n spins- 1/2, its Hilbert space has the ten- 
sor product structure H = Hf n , where Hi = span{|0), |1)} 
is the Hilbert space of a single spin. In the usual computa- 
tional basis we define a reference basis vector 

|0> = |0}i ® ... ® |0)„ (1) 

i.e., all spins up. Let = TVf ™ be the Abelian group of all 
spin flips, where N\ = {1, cr x } acts on a single spin. Ob- 
viously dim H = | A | = 2™ and any vector of the compu- 
tational basis can be written as |i) = .g|0), for some g G N 
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and i = . . . 2" — 1 (as a binary expansion). Moreover, all 
elements satisfy g 2 = 1, g G AT. A generic state in the Hilbert 
space can be written as the superposition of all the possible 
spin flips on this system, namely 

\i>) = a(0)9\O) (2) 

g£N 

witha(.g)eC, E 3 I«(.9)| 2 = 1- 

Consider now the states that can be written as a super- 
position of elements obtained only from a subgroup of spin 
flips G C N acting on |0). Then the corresponding density 
matrix is 

p = \ff,) y,\ = Yl <9)a{g')g\Q)W 

g,g'eGCN 

= a (9ng~g)g\o){o\99 0) 

g,g£GCN 

using the substitution g' = gg in the last equation. 

We now compute the reduced density matrix of an arbitrary 
subsystem A of spins, hence tracing out over all the spins in 
the complement subsystem B. Any element g G G has a 
tensor product structure g = gA® gs, with gA,B £ N acting 
only on A and respectively B subsystems. It is important to 
note that in general gA,B & G. Writing |0) = lO^Os), we 
obtain 

pa = Y a (g)a(gg)gA\OA)(OA\gAgA x 

grgeGQN 

x (0s|gsgsgs|0s) (4) 

where g = c/a® 9b- 

We introduce now two subgroups of G acting trivially on 
the subsystems A and respectively B: 

G A = {geG \ g = g A ®H B } (5) 
G B = {geG \ g = l A ®g B } (6) 

We denote their order by <1a,b = \Ga,b\- With these nota- 
tions, the only non-zero elements in (0B\gBgBgB\0B) satisfy 
9b = 1b (since g 2 B = 11b) and this implies that g G G a- We 
finally obtain 

pa= Y a (gMgg)gA\OA){o A \gA9A (7) 

geG,g£G A 

In general pa will contain off-diagonal terms. The fol- 
lowing lemma gives the necessary and sufficient conditions 
to have a diagonal pa (in the computational basis). 

Lemma 1. The following statements are equivalent: 

(a) pa is diagonal; 

(b) no g G G, g ^ U acts trivially on B, i.e., Ga = {!}; 

(c) no element g in G can be decomposed as the product 
g = 5i ' ff2 with both 31,32 nontrivial and in Ga,Gb 
respectively. 

Proof, (a) <^4> (b): from equation Q Pa is diagonal iff 
cja = 1a, hence Ga — {11} and d,A = 1. 



(b) => (c): Since Ga contains only the identity, there is no 
g% = gA ® H_b different from the identity and thus there is no 
g = gi ■ g2 with nontrivial g\p G Ga.b, which proves the 
sufficient condition. We prove the necessary condition (b) 4= 
(c) ex absurdo. If there were a nontrivial g G G such that 
g = 9 a® 1b, then we can write g — g\- 32 with g\ = g G Ga 
and 32 = 1 £ Gb, contradicting the hypothesis. □ 



B. Entropy of entanglement for a stabilized space 

We are interested to quantify the entanglement present in 
our spin system. Although there is no known entanglement 
measure for a general multi-qubit system, we can study bi- 
partite entanglement of a system described by a pure density 
matrix pab- In this case the von Neumann entropy S is the 
unique measure of bipartite entanglement: 

S = -Tr(pA log 2 pa) (8) 

where pa = Ti b{pab) is the reduced density matrix of the 
sub-system A. The von Neumann entropy of a density matrix 
p is bounded by < S < log 2 d, where d is the dimension of 
the Hilbert space of p. The bound is saturated iff p — l/d, i.e., 
the system is in the totally mixed state. For a bipartite system 
(A, B) of n spins we can readily obtain a simple bound for 
the entropy (using the symmetry S — —Txb{pa log 2 Pa) = 
-Tr A (p B log 2 p B )Y 

< S < min(n J 4, n — tla) (9) 

where tla, n — tla are the number of spins in the A and B 
partition, respectively. 

We now apply the formalism of Section Hi Al to states of 
a stabilized space. Let {U s } be a set of mutually commut- 
ing operators, called stabilizer operators. A state \ip) G H is 
stabilized if is invariant under the action of the stabilizer op- 
erators: U s \ip) = \ip), Vs. Let G be the group generated by 
the stabilizer operators U s . The space £ stabilized by G is 

£ = span{^): U.ty) = \4>)}. 

Suppose now the system is in a state which is an equal su- 
perposition of all elements g G G acting on |0), i.e., a(g) = 
I C I 1 / 2 for all g. This is obviously a stabilized state because 

9\4>) = \Gr 1/2 J2gie G gg'\V) = V 5 6 G. However, 
there are states in the stabilized space which are not an equal 
superposition of all elements g G G acting on |0). Any super- 
position of the form 

\g\-^ 2 Y «WM0) (10) 

heG',g£G 

where G' is a subset of N, is still stabilized by G. 

From now on we will focus on states in the stabilized space 
that are an equal superposition of the elements of G acting 
on the reference state |0). The reduced density matrix for an 
equal superposition state is 

PA = \G\- 1 gA\0A)(0 A \gA~gA (ID 

g£G,g£G A 
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In this case we obtain an analytical formula for the entropy S 
depending only on the boundary of the partition (A, B). 
Define now the quotient G/Gb and let 



f 



_|GJ_ 

\G R \ 



\G\ 
d B 



(12) 



be its order. Notice that / is the number of elements in G 
that act freely on A. If there are I independent generators of 
G acting on A, it turns out that / = 2 l . Define the group 
Gab = G/(G A ■ G B ). We have \G AB \ = \G\/d A d B . 

In the remaining of this section we generalize to an arbitrary 
group of spin flips the results presented in Hill . We can prove 
the following result. 

Theorem 1. Consider a partition (A, B) of the spin system, 
and suppose the system is in an equal superposition of all the 
group elements g £ G C N, acting on the reference state |0). 
The entropy of entanglement is S = log 2 (//dyi) = log 2 — 
\og 2 (d A d B ) = log 2 | Gab | • 

Proof. We first compute the reduced density matrix pa- 
Consider two elements g — gA® g B and g' = g' A <® g' B in G. 
Then g' A — gA if and only if g' = hg, with h £ Gb, and since 
f^ 1 = ds/\G\, from equation il It we obtain 




FIG. 1: A system of spins on an 2D irregular lattice; a typical pla- 
quette and star are denoted by p and s, respectively. The subsystem 
A contains all the spins within the thick boundary. The operators B v 
for the plaquettes situated inside the boundary act only on the spins 
of the subsystem A. The B p of the outside plaquettes that share a link 
with the boundary act on both subsystems A and B and there is one 
such plaquette operator for each link in the boundary of A. Hence 
the entropy is S — Lg A — n c (n c is the number of constraints, see 
text). 



PA = f~ 



E 

g&G A 



9a\0a)(0a\9a9a 



(13) 



Let us compute the square of the reduced density matrix: 



Pa 



= r 2 E 9a\0a}{0a\9a9a9a\0a){0 A \9a9a 



g,g'&G/G B 
g,g'eG A 



f 2 E 9a\0a)(0a\9a9a9 a 



geG/Gs 
g,g'eG A 



= I 2 d A 9A\0 A )(0A\9AgA = f 1 d A PA (14) 

g&G/G B 
§£G A 

Expanding the logarithm in Taylor series we obtain log 2 pa = 
PAfd^ 1 log 2 (dyi//). Then the entropy of entanglement is 



S = log 2 ( fd 



log 2 



Igj 

d A d B 



log 2 \G AB \ (15) 



concluding the proof. □ 
Notice that ifG — N, then N — N A ■ N B and the entropy 
is zero as expected, since in this case the state is an equal 
superposition of all the basis vectors in the Hilbert space. 

Equation dl5> generalizes the result of Ref. II ill (obtained 
for the group of star operators in the Kitaev's model 1 13]) to an 
arbitrary group G of spin flips. We can interpret equation Jl 5i 
as follows. The state of a spin system contains some informa- 
tion. If we have a bipartition (A, B) of the system, we can 
consider the information contained exclusively in A and B as 
the information contained in the bulk of the two subsystems. 
If the system is in a state which is an equal superposition of 
the elements of a (stabilizer) group G of spin flips acting on 



|0), the bulk information is contained in Ga and Gb- The or- 
der of this groups amounts for the "disorder" in the bulk of the 
two subsystems. Then equation dl5t states that the entropy of 
entanglement S is given by the difference between the total 
disorder and the disorder in the bulk. 

Similar results for the entropy of stabilizer states (i.e., for a 
one-dimensional stabilized space), have been obtained in (IF 



C. Spins on a lattice and the entropic boundary law 

The spin systems hitherto considered have no geometrical 
structure and hence we have no notion of what the boundary 
of the two subsystems is. By giving a lattice structure to the 
system, we can find a geometrical interpretation of the equa- 
tion d!5l > as the number of degrees of freedom living on the 
boundary between the subsystems A and B. 

Consider an r-dimensional lattice with n links and a spin- 
1/2 attached to each link (the lattice does not need to be reg- 
ular). Let n p be the number of plaquettes in the lattice. The 
Hilbert space is as before TL = 7if n . Define the stabilizer 
operators 



jedp 



(16) 



acting on the spins belonging to the boundary of any plaquette 
p. The stabilized space is 



£ = m £H: U p \i/>) = Vp} 



(17) 



Let G be the group generated by U p , hence gC = C for 
any g £ G. If all U p are independent, then G is generated by 
the set of all the n. p stabilizer operators. If there are n c con- 
straints on the set {U p }, then the minimal subset generating G 



contains n 



p 



. elements and the order of G is \G\ = 2 np ~ 
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Consider now the stabilized state 

IO = |G|- 1/2 5>|0) (18) 

From Theorem 1, the entropy of entanglement for a state 
|£) corresponding to a partition (A, B) of the lattice is S = 
^(/^a 1 ) = \°g2(\G\/ d Ad B ) = log 2 \Gab\- For this sys- 
tem log 2 6mb) i s m e number of plaquette operators U p acting 
exclusively on A(B). Then S is the number uab of plaquettes 
acting on both the subsystems A and B. 

Hence we can give a geometrical interpretation of equa- 
tion d!5l >. For an r-dimensional lattice, the entropy is equal 
to number uab of degrees of freedom living on the boundary 
between the two subsystems A, B. 

m. THE KITAEV'S MODEL 
A. General formalism 

So far the stabilized states are just some states in a Hilbert 
space and they do not have any physical meaning. Now we 
will analyze a case in which the stabilized states are vectors 
in the ground state manifold of a particular lattice model con- 
structed by Kitaev 11311 . This is a 2-dimensional exact solvable 
spin system on a lattice. Its relevance stems from the fact that 
it was the first example of the new subject of topological quan- 
tum c omp utation II 2lll 311 and because it features topological 
order 0. 

Consider a system of n spins on a (irregular) lattice on a 
Riemann surface of genus g. Again, the Hilbert space is 7i = 
Wf n anddimH = 2". 

The stabilizer operators are the plaquettes 

B P = J] o?, (19) 

jedp 

(j labels all the spins belonging to the boundary of a plaquette 
p) and the stars 

A S =Y[<^ (20) 

jEs 

where j labels all the spins sharing a common vertex s (see 
FigE- On a Riemann surface of genus g the number of sites, 
links (spins) and plaquettes (n s , n and n p , respectively) obey 
the Euler's formula: n s — n + n p = 2(1 — g). By imposing 
also n s = n p , it follows n s = n p = n/2 + 1 — g. We have 
the following two constraints on the stars and plaquettes: 

l[A s = t = l[B p (21) 

Vs Vp 

so there are only n p — 1, n 3 — 1 independent plaquettes and 
stars. 

Let G be the group generated by the n — 2g independent sta- 
bilizer operators {A s ,B p }. We define the protected subspace: 

£ = MeH, A.\tf>) =B p \^) = It/;}} (22) 



so the states in this set are stabilized by G. 
The Hamiltonian of the model is: 

h = -j2 a *-J2 b p (23) 

s p 

The model is exactly solvable because all the stabilizer oper- 
ators commute with each other (since they share either or 2 
links) 

[A s ,B p }=0, Vs,p (24) 

Its ground state is the protected subspace manifold C. 

We now show that the ground state manifold L for a genus- 
g Riemann surface is 2 2s -fold degenerate. Since all the stars 
and plaquettes commute, we can label n s + n p — 2 = n — 2g 
spins out of n and thus the dimension of the ground state is 

dim/: = 2"- ( "- 2b) = 2 2fl (25) 

hence the system exhibits topological order fl4ll . 

Another way to see the same thing is to notice that this 
model features string condensation |14]. Let j z (j x ) be a 
curve connecting sites along the links of the lattice (dual lat- 
tice) as in Fig.|2(for a square lattice). We can define two types 
of string operators (simply called "strings"): 

(i) a z-string is the product of all a z operators along the links 
belonging to a curve j z running on the lattice; 

(ii) an x-string is the product of all a x operators along the 
links crossed by a curve ^ x (hence running on the dual lat- 
tice). The action of an x-string is to flip all the spins (i.e., 
links) intersected by the curve r y x . 

More formally, a string operator is: 

= n °i > a = x > z < 26) 

j'67 a 

By j £ 7 Z (j G j x ) we mean all the links belonging to 
(crossed by) the string -f z (j x ). A string-net is a product of 
string operators. 

Closed strings of both types commute with the Hamiltonian 
O, 

[W a [-t a ],H]^0, a = x,z (27) 

where 7° is a loop on the lattice or on the dual lattice. This 
is because a closed string has either or 2 links in common 
with any plaquette or star. However, open strings do not com- 
mute with the Hamiltonian. More precisely, open z-strings 
(x-strings) anticommute, with star (plaquette) operators with 
which they share (cross) a single link (see Fig.|5J- 

Equation i27\ implies that there are closed strings in the 
ground state of the Hamiltonian d23l >. We say that we have 
closed string condensation 1 14] in the sense that closed strings 
are present in the ground state and that they are not made of 
smaller pieces that are still present in the ground state. Smaller 
pieces are indeed open strings and they pay energy. The group 
G generated by the closed strings is the group of all the closed 
string-nets of x- and z-type. 

On a genus-g Riemann surface (g > 0) there are con- 
tractible and non-contractible loops. A loop on the lattice is 
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FIG. 2: (Color online) A k x k square lattice of the torus; opposite 
boundaries are identified. The end of an open z-string anticommutes 
with the star based at the site where the string is. Similarly, the end 
of an open x-string anticommutes with the plaquette on the square 
where it ends. The figure also shows that a star on the lattice corre- 
sponds to a plaquette in the dual lattice and vice versa. 

contractible if is homotopic to the boundary of a plaquette. 
It turns out that all the string operators based on contractible 
loops are made of products of stars and plaquettes. This also 
implies that string operators based on contractible loops have 
a trivial effect on the ground state |£) e C: 

tnr]l£) = I]>;IO = IO (28) 

ier 

where T is a contractible loop on the lattice (or the dual lat- 
tice). 

Consider now the non-contractible loops. It is enough to 
consider only string operators associated to non-contractible 
loops with winding number 1, [Since all string operators 
square to the identity, string operators associated to non- 
contractible loops with winding number n are equal (modulo 
a product of stars) to the ones with winding number n mod 2]. 
The associated string operators cannot be written in terms 
of products of star and plaquette operators and hence they 
have a nontrivial action on the states. But since they still 
commute with the Hamiltonian, they map ground states into 
ground states. The algebra L(£) of linear operators acting 
on the ground state manifold C is the algebra of the closed 
string operators of x- and z-type. However, the contractible 
string operators have a trivial effect on £, so only the non 
contractible ones matter. Consider the string operators as- 
sociated to the non-contractible loops {7;,?' = 1, ...,4g} = 
{7i> -!72 B . 7! B +i ! --:71b}; the loops { 7 f , ...,7f B } generate 
the homotopy group of the Riemann surface. We label the 
loops such that 7? and 7 4 ? + 9 fl intersect, with i = 1, .., 2g, see 
FigEl 

Wi = W a bi] i = l,...,4 fl (29) 

where a — x for i = 1, 2g and z otherwise. We refer 
to Wi,i = 1, ...,2g as "ladder" operators since they flip all 



the spins along a ladder going around the non contractible 
loops of the surface. The pair (w.i, w;+2g) has the same com- 
mutation relations as (<r x , a z ) and generates a 4-dimensional 
algebra. We see that for any i = 1, 2g we have a copy 
of the same algebra. Then we have 2g mutually commuting 
copies of the same algebra a x ,a z and hence L(£) is 4 2fl - 
dimensional: 

w t w i+2s = -w i+2g Wi, i = 1, 2g 
[wi,Wj]=0, j^i±2g, i = 1,...,4q (30) 

Therefore the ground state manifold C is 2 2fl -fold degenerate. 
This degeneracy is the sign of the topological order of this 
system Il4ll . and is robust against arbitrary perturbations fl3l 
Y\M . Topological order is the notion needed to describe those 
states of the matter like fractional quantum Hall liquids fl7ll 
which are not explained by the Landau theory of symmetry 
breaking with local order parameters 0]. 

We want now to give an explicit expression for the states in 
the ground state manifold. Let A be the group generated by 
the stars A s . Let T be the group generated by the 2g ladder 
operators of x-type Wi, with i = 1, 2g. Then the elements 
of T are of the form 

w(s) = Y[w s / (31) 

j=i 

where s = (si,..,S2 g ) and Sj = 0,1, which implies that 
\T\ = 2 2s . We will call N the group of all closed string-nets 
of the a;- type: 

N = A • T (32) 

For a generic lattice with n s stars, n spins and n p plaque- 
ttes, the number of independent star operators is n s — 1 = 
n — n p + 1 — 2g and hence the order of A is \A\ — 2 Tls ~ 1 = 
2 n ~ n p+ 1 - 2 s. The order of N is \N\ = 2™-- 1+2 s = 2™"™f +1 . 

Since on the lattice any loop intersects the boundary of a 
plaquette in an even number of points, it follows immediately 
that the states u;(s)|0) are stabilized by the plaquettes B p . 
Then we have 2 2fl states in the stabilized subspace C given 
by 

|£(s)> = |,4|-Vau,(s) £ g\0) = w(sM(0)) (33) 

g&A 

which are mutually orthogonal by construction. This shows 
again that the ground state manifold is 2 2b degenerate. Each 
of these states is an equal superposition of the elements in A 
and so it falls under the hypothesis of Theorem 1 . Notice that 
an arbitrary superposition of the |£(s))'s is still a ground state, 
but obviously it is not an equal superposition of the elements 
of a group: 

10= E a ( s M s )2l°> (34) 

g&A,w(s)eT 

so C = span{|£(s)}} = span{T|£(0))}. 
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FIG. 3: (Color online) A lattice region consisting of all the spins in- 
side or crossed by a loop (thick line) on the dual lattice. The Ea 
stars based on the sites inside the region A (red diamonds) act ex- 
clusively on the subsystem A and the E_g stars based on the black 
sites act only on the subsystem B. There are Eab stars (white dots) 
acting on both subsystems A and B. The total number of stars is 
n a — Ea + Es + Eab = n — n p + 2 — 2g, where n (%,) is the 
number of spins (plaquettes) of the lattice. 



We now prove that all the basis states |£(s)) have the same 
entanglement. 

Proposition 1. For a given lattice partition (A, B) all the 
states |£(s)) have the same entropy of entanglement S. 

Proof. We can decompose the ladder operators as w(s) = 
wa(s) ®wb(s), where Wa(b){s) acts only on the A (B) sub- 
system. From the circular property of the trace and w 2 (s) = 
H, w'(s) — w(s), s = 0, 2 2fl — 1, it follows immediately 
that all the basis states |£(s)) have isospectral reduced den- 
sity matrices. Since |£(s)) = io(s)|£(0)), then /9a[£( s )] = 
Tr B (|£(s))(£(s)|) = Tr B Ms)|£(0))<£(0)Hs)) = 
Wa (s)Trg (|£(0)) (£(0) \)wa (s) . Therefore 

S(|€(s)» = S(|£(0)», Vs = 0...2 2 *-1 (35) 

and all basis states have the same entropy. □ 
Let us now compute the entropy of entanglement for a state 
|£(s)). The number of closed string-nets acting exclusively 
on A(B) is <1a.b — 2 Ea b , where T>a,b is the number of in- 
dependent star operators acting exclusively on the subsystem 
A(B). From Theorem 1 we obtain 

2n s — 1 

S = lo S2 , , = Us ~ 1 - - S S 
dACLB 

= n-n p + l-2g-Z A -X B (36) 

We notice that the topological order in this model manifests it- 
self in both the degeneracy and the entanglement in the ground 
state. This suggests the very appealing possibility that entan- 
glement could detect topological order. 

The total number of lattice sites is n s = + £ B + ^ab = 
n — n p + 2 — 2g, where 'Sab is the number of independent 
star operators acting on both subsystems A and B. We obtain 

S = E AB - 1 (37) 



If we choose the partitions in a convenient way, we can give 
a clear geometrical picture of the formula (I37> . Let A be the 
set of all spins inside or crossed by a contractible loop in the 
dual lattice (see Fig|5}. The spins intersected by the loop are 
the boundary of A, while the ones inside are the bulk. If the 
loop is convex, the number of spins ul in the perimeter L of 
A is til = Sab (see Fig|3jl. Therefore 

S = n L -l (38) 

It is interesting that no partition has zero entanglement for 
all the |£(s)) states. The argument is simple: S = <^> \N\ — 
dAds A = Aa ■ Ab', but this cannot be satisfied, since 
there is at least a star or a ladder acting on both A and B for 
any partition (A, B), hence S > 0. The group A splits in 
A = A A ■ Ab for any partition only if is the group generated 
by the single spin flips, namely N. If every spin is shared by 
at least two generators of A (which is always the case for star 
operators on a lattice), then the entropy cannot be zero and we 
have an absolute entropy. 



IV. GROUND STATE ENTANGLEMENT FOR THE 
KITAEV'S MODEL ON A TORUS SQUARE LATTICE 

In this section we consider a square k x k lattice on the torus 
(0 = 1) and we calculate explicitly the entropy S for several 
bipartitions (A, B) of the lattice. On such a lattice there are 
n s = k 2 sites and n — 2k 2 spins. 

We have two ladder operators W\ and w-2 corresponding to 
the two non-contractible loops which run along the parallel, 
and respectively the meridian, of the torus. The group A gen- 
erated by the stars has order \ A\ — 2 n > 2 ~ 1 . Then the group TV 
is generated by A, Wi and ^2 and hence |iV| = 2"/ 2+1 . 

The ground state is four-fold degenerate and the vectors 
ii 3 = 0) 1 form a basis, with 

|&> = M- 1 ^ t4 5>l°> (39) 

g€A 

An arbitrary vector of the ground state |£) 6 C can then be 
written as 

l 

If) = H a ij\&j) 
i,3=0 

= \A\~ 1/2 ^(aoo + aoiwi + a w w 2 + a 11 w 1 w 2 )g\0) 

g&A 

= lAl-^^UglO) (40) 

g&A 

where U = aool + aoiWi + cliqW2 + auWiW2 and 

The associated density matrix is 
l 

P= aijaimwlw^pow^wl (41) 

i,j,l,m=0 
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Any element g G A leaves invariant the ground state, 
ffl&j) = l&i). hence 5 A) = Po. where 

p = ICooXCoo (42) 

There is another important issue to point out. Suppose we 
have two ladder operators wi and w\ with homotopic sup- 
ports. This means that they are related by an element g 6 A, 
W\ = gw\. Since gp = p a and [g, wi] = 0, then 

wipo = gwip = wipo (43) 

Therefore both w\ and wi have the same effect on po and 
hence we can work only with a representative w\ . 

A proposition useful for computing the entropy in some of 
the examples below is the specialization of the diagonality 
condition for the general Kitaev model. We first prove that 
a generic product of spin flips acting exclusively on the sub- 
system A commutes with the Hamiltonian only if is a product 
of closed strings. Let N a be the set of closed string nets that 
act exclusively on A, N a = {g G N : g = gA <8> Is}. Then 
we have the following: 

Lemma 2. Let gA <& Is G Na be a generic product of spin 
flips operators acting exclusively on A. Then [H, gA <S> is] = 
<£> ® is G ATa. 

Proof. The proof is obvious since only closed string nets 
commute with the Hamiltonian, so if gA ® Is commutes with 
H it must be in N and henceforth in iV^. □ 

Proposition 2. Suppose the system is in the ground state 
£oo) given by equation J39I . Then the reduced density matrix 
Pa is diagonal (i.e., cIa — 1) if and only if for every g G 
A, gA 7^ 1U, 9a does not commute with the Hamiltonian: 
[9A,H}^0. 

Proof. "=>•": Let d,A — 1, and suppose, ex absurdo, that 
there is a non trivial element gA G N such that [gA, H] = 0. 
Then from Lemma 2 this means that gA <S> Us is a nontrivial 
element of .Aa, contradicting the hypothesis. 

The reverse implication "-*=" obviously holds because for 
any nontrivial gA which does not commute with the Hamil- 
tonian, then it follows from Lemma 2 that g = gA ® Us is 
not in Aa, and thus only the identity belongs to Aa, which 
concludes the proof. □ 

A corollary is that if the system is in the ground state 
\Cij), the same proposition holds true by substituting A with 
wlw^A. If the system is in a generic ground state one just 
replace A with N. 

In the following we compute the entropy for several sub- 
systems A for both |£oo) and the generic |£) ground states. 
Although the general ground state |£) is no longer an equal 
superposition state and hence we cannot apply Theorem 1, in 
some cases we can calculate explicitly the von Neumann en- 
tropy for For completeness we also review some examples 
given in Ref. full for the ground state |£oo)- 

A. One spin 

Take A to be a single spin. In this case it is obvious that 
no closed string net g £ N acts exclusively on A, hence pa 



(a) (b) (c) (d) 

FIG. 4: (Color online) The subsystems A (thick/blue spins) used in 
calculating the entropy S: (a) the spin chain; (b) the vertical ladder; 
(c) the cross; (d) all the vertical spins. 



is diagonal (from Lemma 1) and both eigenvalues are equal to 
1/2 (from symmetry, the entries for spin-up and spin-down are 
equal). Since the entropy is S = 1, it follows that any spin is 
maximally entangled with the rest of the system. 

B. Two spins 

We want now to compute the entanglement between two 
arbitrary spins of the lattice. To do this, we first obtain the 
reduced density matrix pij of the two spins by tracing out all 
the other spins. Since we want to calculate the entanglement 
between the two spins, we use as an entanglement measure the 
concurrence C of the mixed state p^ of the two qubits (i.e., 
the two spins) defined as fl9ll : 

C = max{0, ^/\[ - ^ - v/As - (44) 

where Ai, A2, A3, A4 are the eigenvalues (in decreasing order) 
of the matrix p t3 {cr v ® a y )p* J (a y ® a v ) and a y is the Pauli 

matrix. Since no closed string net g G N acts exclusively 
on the two spins, the reduced density matrix is diagonal 
(again from Lemma 1). Let = diag(a, b, c, d). A simple 
calculation shows that C = always, hence there is no two- 
qubit entanglement between any two spins. 

We see that although an arbitrary spin is maximally entan- 
gled with the rest of the system, the entanglement is zero be- 
tween any pair of spins. 

C. The spin chain 

Let A be the set of the k spins belonging to the meridian 
7 2 i of the torus as in Fig|4ja) and consider the system in the 
state |£oo)- Since the state |£oo) is the equal superposition of 
all the group elements in A acting on |0), we can apply The- 
orem 1 and the entropy is S = k 2 — 1 — log 2 dA — log 2 c?s- 
It is obvious that any spin flip on the chain does not com- 
mute with the plaquettes sharing that spin, so from the Corol- 
lary to Proposition 2 no closed string-net g G N acts exclu- 
sively on A and thus pa is always diagonal for any ground 
state In particular, for |£no) we have dA — L Therefore 
S = log 2 (|G|/ds) = log 2 / from equation H2\ . The number 
of possible configurations of spins on the chain 7 2l is 2 k , but 
there are only / = 2 k ~ 1 different configurations of spins in 
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A that enter the ground state, namely the ones with an even 
number of spin flips. Indeed, we have k — 1 stars acting in- 
dependently on the chain and we can obtain all the allowed 
configurations applying products of these stars (i.e., elements 
in A acting freely on the chain), which gives / = 2 k ~ 1 . Then 
the entropy is 

S = log 2 f = k-l (45) 

Let now the system be in a generic |£) of the ground state 
manifold. This state is a superposition with arbitrary coeffi- 
cients of the four orthogonal states in C: 

l l 

io = £ atfi&> = \a\-^ 2 J2 <*i<*>\ E ^ 

i,j—0 i,j=0 g£A 

Since this is not an equal superposition of elements of a group 
(acting on |0)), we cannot apply Theorem 1. Nonetheless, as 
shown before, the reduced density matrix for this system is 
diagonal (in the computational basis). Moreover, all the pos- 
sible configurations on the chain are allowed in A^O) since 
all of them can be realized by applying some horizontal lad- 
der w. Thus there are 2/ = 2 k configurations of the chain 
that are in N\0), where / = 2 k ~ 1 is the number of even 
(odd) spin flips configurations. We see that the states |£oo) 
and |£iq) give states with an even number of spin flips on the 
chain, while |£oi) and |£n) give states with an odd number 
of spin flips. The eigenvalues corresponding to the even spin 

flip configurations are then |^l|<iB(|aoo| 2 + |aio| 2 ) = |^4|rfsQ!, 
while the eigenvalues corresponding to an odd number of spin 
flips on the chain are |^|ds(|a i| 2 + |an| 2 ) = |^4|rfs(l — 
a). If \i are the previous eigenvalues, the entropy is S = 

~ Yh=i A. lo g2 Ai = fc - 1 - alog 2 a - (1 - a) log 2 (l - a). 
Defining the binary entropy 

H 2 (x) = -x\og 2 x-(l-x) log 2 (l - x) (47) 

we obtain 

S = k - 1 + H 2 (a) (48) 
D. The spin ladder 



since the ^4-system is in the totally mixed state. 

What happens if the system is in a generic ground state 
In general, the reduced density matrix is no longer diagonal 
(in the computational basis). From equation J4-1I we find 

l 

Pa = E a ija; m Tr B (^u> 2 pow™ui 2 ) (51) 

Since the set A is the vertical ladder, both ladder operators 
have a particularly simple action: w\ (w 2 ) acts only on sub- 
system B (A). Then Tr B (w 2 po) = W2Tib(po)- Moreover, 

Tr B (wip ) = J2g.g'eA X A\ Q A)(0A\xAX l A (0B\x' B W 1 \0 B ). 

These scalar products are different from zero if and only if 
x' B = wi\. This would imply that x' A ®x' B — x' A ®W\ isacon- 
tractible string net in A and this can happen only if x' A is a lad- 
der operator acting fully on A, which is impossible (notice that 
a double ladder is a product of stars and hence a contractible 
string net). Thus Tvb{wiPq) = 0. Similarly, Txb{wiW2Pq) = 
w 2 Tr B {wipo) = 0. We also have Tr B {wip wi) = Tv B (p ) 
and Tr B {w2Pow-2) = W2^b{po)w2- From equation ( l50l we 
know that for the ladder Tib(pq) = 2~ k i A and we obtain 

p A = 2- k {t + pw 2 ) (52) 

with p — 2Re(aooaio + aoian)- Since iwf = 1 and 
Tr(w 2 ) = 0, the eigenvalues of w 2 are ±1 and they have the 
same multiplicity, namely 2 k ~ 1 . Hence the eigenvalues of pa 
are A± = 2~ fc (l ±p) and the entropy is 

S = k-1 + H 2 (— ^) (53) 



E. The cross 

The subsystem A includes all the thick (blue) spins in the 
state FigHc). This is a system of 2k spins. Let the system be 
in the state |£oo)- Again, no element of A is able to flip spins 
only on this subsystem, so <Ia = 1 and the reduced density 
matrix pa is diagonal (in the computational basis). There are 
k 2 — 1 — (2k — 1) stars acting independently on B. The entropy 
is thus S = 2k - 1. 



Let A be the set of spins belonging to a "ladder", i.e., all 
the spins intersected by the curve 7^1, Fig|4jb). Again, the 
A partition contains k spins. Consider the system in the state 
|£oo)- From Proposition 2, no g G A acts exclusively on the 
subsystem A, hence d A = L The number of independent stars 
acting only on B is k 2 — 1 — 2k + k (there are k 2 — 1 — 2k 
stars which do not touch the subsystem A, plus k pairs of stars 
based on the two ends of each spin in A, leaving it invariant). 
Then d B = 2 k ~~ k ~ 1 and the entropy is 

S = k (49) 

This implies immediately that 

Tr s (po) = 2~ k l A (50) 



F. The vertical spins 

We now take A to be the set of all vertical spins of the lat- 
tice; then B is the set of all horizontal spins, see Fig|4jd). The 
system is considered in the state |£oo)- Since in this case no 
closed string operator g 6 A acts trivially on either subsys- 
tem, we have dA = d B = 1. The entanglement is S = k 2 — 1, 
which is the maximum possible value for a state. 

G. The disk 

In this section we take the system A to be a disk, i.e., a 
region homeomorphic to a plaquette. Let us start, as usual, 
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(a) (b) 



FIG. 5: (Color online) A region A of the lattice obtained by taking 
all the spins (thin/black lines) inside or crossed by a loop (thick/red 
line). £a,b are the number of sites whose stars act exclusively 
on A (diamonds/red sites) and respectively, B (solid/blue circles); 
<1a,b = 2 Ea b . The number of sites with stars acting on both 
subsystems (open/white sites) is T,ab = ni + 712 + 713; m is the 
number of (open/white) sites having % nearest neighbours inside A, 
Area conservation implies Ea + Es + Eab = k 2 . The entropy is 
S = Eab — 1. (a) If the boundary is a convex loop (i.e., a rectangle), 
Eab is equal to the boundary length Lqa (in lattice units), since in 
this case 712 = n 3 =0; the entropy is S — Lqa — 1. (b) If the 
boundary is an arbitrary loop, Lqa = ni + 2ri2 + 3ri3, hence the en- 
tropy is S = Lqa — 712 — 2ri3 — 1; in the figure 712 = 1 (open/white 
triangle) and n 3 = 1 (open/white square). 

by assuming the system is in the state |£oo}- Consider a con- 
tractible loop 7 on the dual lattice and let A be the set of spins 
inside or intersected by 7 (Fig. [5}- Let be the number of 
sites whose stars act only on A, respectively B; Y±a,b is the 
area of A, B in lattice units. Let Sab be the number of sites 
with stars acting on both subsystems. We also have d A .B — 
2 Sa < b . Area conservation implies + + Y>ab = k 2 . 
The entropy is 5" = k 2 - 1 - S A - E B = T, AB - 1. Let 
ni, i = 1...3, be the number of sites in Y>ab having i nearest 
neighbours inside A. Then 

^ab = n x + n 2 + n 3 (54) 

and the boundary length is 

Loa = rii + 2n 2 + 3rt 3 (55) 

If 7 is a convex loop (Fig. |5^), then Y>ab = Lqa (since 
fii = ns — 0) and the entropy is equal to the perimeter of the 
boundary (up to a constant) S — Lqa — 1. In general, if the 
boundary of A is an arbitrary loop 7 on the lattice (Fig.|5p), 
the entropy is 

S = L aA - n 2 - 2n 3 - 1 = ^( l qa + 2n x + n 2 ) - 1 (56) 

We see that for the disk S is exactly the geometric entropy 
l2(ill of a spatial region A. 

It is easy to prove that the entropy is bounded from above 
and below by two linear functions. From eq. d55l > it follows 
that rii < Lqa/i and 2n 2 +3n3 < Lqa- Then eg. d56l implies 

\lq A - 1 < S < 7 -L dA - 1 (57) 



TABLE I: The entropy S of the systems analyzed in text for two 
ground states, |£oo) an d the generic for two spins the value 
shown is the concurrence C. The constants are a — \aoo\ 2 + \aio\ 2 
and p = 2Re(aooSio + ffloian); H 2 (x) = -a;log 2 a; — (1 — 
x) log 2 (l — x) is the binary entropy. 
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two spins i,j 


c = o 


c* = 


2. 


spin chain 


k- 1 


k-l+H 2 (a) 


3. 


spin ladder 


k 


k-l + H 2 (^) 


4. 


cross 


2k - 1 




5. 


vertical spins 


k 2 -l 




6. 


the disk 


LaA -TL2- 2n 3 - 1 


LaA -TI2- 2n 3 - 1 



Let us now consider the system in the generic ground state 
£}. The expansion of the reduced density matrix is pa = 
i,m=o ^i3^im^B{w{wlpQwf'w 2 ). From equation (03j, 
for the disk we can choose the ladders wx,2 such that they 
act only on B, hence WiA — ^a< * — L 2. Then it fol- 
lows immediately that Tr B(wiPoWi) = Ttb(po),i = 1,2, 
and Ttb(vjiw 2 pqWiw 2 ) = Tib(po)- As in section llV Dl 
Ti"b(po w i) = 0, since the ladder W\ cannot act exclusively on 
subsystem A. A similar reasoning implies also Tib(pow 2 ) = 
= Tib(poWiw 2 ) and we obtain pa = Tib(po)\ hence the 
entropy is the same as in the previous case, S = Lqa — n 2 — 
2ns — 1. Thus for the disk the entropy obeys the boundary law 
for any ground state 

V. CONCLUSIONS 

Apart from being one of the most striking conceptual fea- 
tures of quantum mechanics, entanglement proves also to be 
a powerful tool in the study of many body spin systems, as 
several articles pointed out recently. 

The first topic to which our article is related is the study 
of entanglement in spin systems. Several authors have calcu- 
lated the entanglement in ID spin chains. In the case of XY 
and Heisenberg models, the authors in Refs. |7, 8] calculated 
the entanglement between a spin block of size L and the rest 
of the chain. They found two characteristic behaviors. For 
critical spin chains, the entanglement scales like S ~ log 2 L, 
whereas for the noncritical case S saturates with the size L of 
the block. This result is in agreement with the result for black- 
hole thermodynamics in 1+1 dimensions i2lll22Tl . which sug- 
gested a connection between the entanglement measured in 
quantum information and the entropy of the vacuum in quan- 
tum field theories. 

In this article we investigated bipartite entanglement in spin 
systems for states in a stabilized space. For a bipartite system 
in a pure state, the von Neumann entropy of the reduced den- 
sity matrix is the unique measure characterizing the entangle- 
ment between the two subsystems. We showed that for states 
that are an equal superposition of all the elements of a stabi- 
lizer group generated by spin flips, the entanglement entropy 
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of a bipartition (A, B) depends only on the degrees of free- 
dom belonging to the boundary between the two subsystems. 
This property provides an interesting link to the holographic 
principle. As an example we studied the entanglement present 
in the ground state of the Kitaev's model. Apart from its spe- 
cial interest in quantum computation (it was the first example 
of topological quantum computing), this model is also rele- 
vant per se, due to the nontrivial topology and to the specific 
nature of the spin-spin interaction which generates topological 
order. On a Riemann surface of genus g the degeneracy of the 
ground state is 4 fl and it is stable against local perturbations. 
We found analytical results for the ground state entropy of 
several bipartitions (A, B) of a toroidal square lattice. In this 
case, although no two spins of the lattice are entangled (the 
concurrence is zero for any pair of spins), the ground state has 
genuine multi-body entanglement. For a convex region A of 
the lattice, its geometric entropy is linear in the length of the 
boundary. Moreover, for states which are an equal superpo- 
sition of all elements g E G C N acting on |0), no partition 
has zero entanglement, so the system has an absolute entan- 
glement. Finally, we argued that entanglement can probe the 
topology of the system and raised the very interesting question 
of whether it could detect (or measure) topological order. 

It is relevant to put our results in perspective and to com- 
pare them with known results. The holographic principle 
(HP) emerged recently as a paradigmatic universal law 11511 . 



A simplified statement of HP is: The maximum entropy of a 
region is proportional to the area of its boundary. It appar- 
ently contradicts the naive expectation that the entropy of a 
region should be proportional to its volume. The entropic area 
law appears as a recurrent pattern and has been recovered in 
several (apparently unrelated) physical systems. In black hole 
thermodynamics it is expressed as the Beckenstein-Hawking 
law, Sbh =A/A: the entropy of a black hole is a quarter of 
its horizon area (in Planck units) 1!23[| . For a scalar field in 
2+1 and 3+1 dimensions, Srednicki 1 24] showed that the en- 
tropy of a region R is proportional to the area of its boundary, 
and not to its volume. Recently Plenio et al. 12511 found an- 
alytically the same behavior for the entropy in the case of a 
harmonic lattice system in d-dimensions. 

The entropic boundary law recovered in this article for spin 
systems provides another instance of the universality of HP. 
Due to the close relationship between concepts like entropy 
and entanglement, the holographic principle gives insights 
into fundamental questions of quantum information theory, 
like What is maximum information content needed to describe 
a region R? or How much information can be stored inside 
a system A ? This confluence of diverse fields, like black hole 
thermodynamics, QIT, and spin systems, can bring together 
insights and shed new light on fundamental problems. 

PZ. gratefully acknowledges funding by European Union 
project TOPQIP (contract IST-2001-39215). 
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